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Abstract 

NASA Glenn Research Center (GRC) and the Politecnico di Milano (POLIMI) have initiated a joint 
propagation campaign within the framework of the Alphasat propagation experiment to characterize rain 
attenuation, scintillation, and gaseous absorption effects of the atmosphere in the 40 GHz band. NASA 
GRC has developed and installed a K/Q-band (20/40 GHz) beacon receiver at the POLIMI campus in 
Milan, Italy, which receives the 20/40 GHz signals broadcast from the Alphasat Aldo Paraboni TDP#5 
beacon payload. The primary goal of these measurements is to develop a physical model to improve 
predictions of communications systems performance within the Q-band. Herein, we describe the design 
and preliminary performance of the NASA propagation terminal, which has been installed and operating 
in Milan since May 2014. The receiver is based upon a validated Fast Fourier Transform (FFT) l/Q digital 
design approach utilized in other operational NASA propagation terminals, but has been modified to 
employ power measurement via a frequency estimation technique and to coherently track and measure 
the amplitude of the 20/40 GHz beacon signals. The system consists of a 1 .2-m K-band and a 0.6-m Q- 
band Cassegrain reflector employing synchronous open-loop tracking to track the inclined orbit of the 
Alphasat satellite. An 8 Hz sampling rate is implemented to characterize scintillation effects, with a 1-Hz 
measurement bandwidth dynamic range of 45 dB. A weather station with an optical disdrometer is also 
installed to characterize rain drop size distribution for correlation with physical based models. 


I. Introduction 

Though the primary focus of NASA’s future space communications architecture is in the exploitation of Ka 
and optical frequencies, NASA is also investigating the potential use of spectrum in the Ka/Q-bands (37- 
42 GHz) and V/W-bands (74-84 GHz) as a downlink option in the next generation Space Based Relay 
(SBR) that is expected to replace the existing Tracking and Data Relay Satellite (TDRS) system in the 
2025 timeframe. As such, NASA Glenn Research Center is spearheading a small, collaborative effort to 
characterize atmospheric propagation effects at these millimeter wave frequencies. A major part of this 
campaign includes the use of the Alphasat Aldo Paraboni TDP#5 20/40 GHz beacon in Europe. The 
primary goal of these measurements is to develop a physical model to improve predictions of atmospheric 
attenuation within the desired spectrum. Herein, we describe the design and performance of the 20/40 
GHz beacon receiver, built in-house by NASA Glenn Research Center, and which has been recording 


data at the Politecnico di Milano (POLIMI) campus since May 2014. The receiver is based upon the 
validated FFT digital design approach utilized in other operational GRC propagation terminals [1-2]. The 
new design incorporates upgrades and modifications to coherently track and measure the amplitude of 
the 20/40 GHz beacon signals with improved dynamic range. It is believed that these measurements, 
combined with thorough characterization of local meteorological conditions, will provide the necessary 
basis for improved models to assess the performance of communications systems employing frequencies 
in the Q-band. 


II. GRC K/Q-band Beacon Receiver Design 

The Alphasat beacon receiver developed at NASA GRC consists of a 1.2m K-band and a 0.6m Q-band 
Cassegrain reflector with equivalent antenna beamwidths of 0.9 deg. The receivers employ independent 
open-loop tracking systems for each antenna to track the inclined orbit of the Alphasat satellite. A block 
diagram of the front ends and common Intermediate Frequency (IF) downconversion cards (labeled K- 
DACS and Q-DACS) of the system is provided in Figure 1 and Figure 2, respectively. Primary 
downconversion to a common IF of 70 MHz is performed within the temperature controlled RF boxes 
mounted directly behind the antennas. Independent temperature control of the LNA and the plate upon 
which the RF electronics are mounted maintain a temperature stability of +/- 0.01 deg C. From the RF 
boxes, the signal is routed to a secondary temperature controlled IF box where final downconversion 
stages take place to the common 455 kHz IF. Temperature stability of the IF box is maintained to within 
+/- 0.25 deg C. 
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Figure 1 . Block diagram of coherent K/Q-band beacon receiver 
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Figure 2. Block diagram of 70 MHz to 455 kHz downconversion card (labeled K-DACS and Q-DACS in Figure 1). 






The IF Cards (K-DACS/Q-DACS) are integrated circuit boards originally utilized during the ACTS 
experiment which were subsequently modified to maintain coherent downconversion of the receive 
beacon signals. The input to these cards necessitate a 70 MHz IF and provide the necessary filtering and 
final downconversion stages to the ultimate 455 kHz IF. To accomplish this, the 19.701 GHz K-band 
beacon is first downconverted in the RF box via a single stage to the 70 MHz IF frequency input. Due to 
the limited availability of components (i.e., local oscillators and filters), the 39.402 GHz signal is 
downconverted in two stages to obtain the necessary 70 MHz IF signal input to the Q-DACS card. All 
downconversion stages are referenced to a common ultra-stable 10 MHz reference oscillator housed 
within the IF box. 

From there, the signals are routed via coaxial cable to the control and data acquisition computer. The 
455 kHz IF signal is Nyquist sampled by a 12-bit National Instruments 5124 data acquisition (DAQ) card 
at a sampling frequency (fs) of 1.111MHz. At an 8 Hz data measurement rate, 2 17 samples are collected 
for a final fs/N resolution of 8.47Hz. The 8 Hz power measurements are averaged each second to derive 
1 Hz data and both data sets are recorded to file. Frequency estimation and power measurement of the 
beacon signal employs a modified Quinn-Fernandes frequency estimation routine, as described in [3] and 
discussed in more detail in the next section. During deep rain fades, i.e., fades greater than 30dB, the 
beacon receiver utilizes the K-band signal for frequency tracking to maintain lock on the Q-band signal. 
Based on this approach, the measurement range of the Q-band receiver is improved down to a 3dB bin 
signal to noise ratio (SNR) and is capable of reestablishing lock once the signal emerges from the noise. 
This implementation provides a measurement dynamic range of approximately 38 dB for the 8 Hz 
sampling rate, and 45 dB for the 1 Hz samples. Further, by employing a frequency estimation routine, we 
are also able to derive Doppler parameters for estimates of relative spacecraft position and calibration of 
gain variation across the passband of our receiver. 

The receivers were installed atop the roof of the Dipartimento di Elettronica, Informazione e Bioingegneria 
(DEIB) building on the POLIMI campus in April 2014 and is shown in the photograph of Figure 3. During 
initial operations, several anomalies were addressed and operational status was achieved on June 1, 
2014. To supplement the beacon receivers, a suite of weather sensors are utilized to obtain surface 
temperature, pressure, humidity, wind speed/direction, and rain rate. A disdrometer has also been 
installed at the POLIMI site to obtain drop size distribution, fall velocity, and high resolution rain rate 
measurements during rain events. 



Figure 3. Photograph of K/Q-band beacon receiver installed at POLIMI campus. 



III. Frequency Estimation Approach 

For a digital FFT receiver in the absence of a phase-locked loop (PLL), the beacon frequency can drift 
and must be identified such that a power measurement can be performed. As the frequency drifts from 
the center of an FFT bin, the signal power will spill over into adjacent bins and, if measuring the peak bin 
power, will decrease by as much as 3 dB at the bin edges. This effect can be reduced by summing over 
multiple FFT bins, but this has the negative effect of reducing the system dynamic range due to the 
integration of the noise power over a wider bandwidth. One method to alleviate this problem is to employ 
a frequency estimator to more accurately identify the signal frequency within the peak bin and measure 
the signal power at baseband. 

Figure 4 shows the simulated performance of various frequency estimators as compared to the single bin 
power estimation technique for an SNR of -10 dB [3], As the simulated frequency is swept between 
adjacent FFT bins, we observe that any of the identified frequency estimator routines sufficiently track the 
frequency and provide an accurate power measurement, as compared to the peak FFT approach, where 
the scalloping effect is readily observable. 
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Figure 4. Frequency estimate of various estimation routines (top) and the corresponding power (bottom) for a simulated -10dB SNR. 


For our receiver design, we employ the Quinn-Fernandes frequency estimation technique [4], The Quinn- 
Fernandes technique attempts to find the local maximum of a smoothed periodogram nearest some initial 
frequency estimate by optimizing the parameters a and p in the filtered spectrum as, 

y t - Py t - 1 + y t - 2 = £ t- 1 + £ t-z 

Utilizing this approach, it is possible to derive the frequency estimate within typically 2-3 iterations, 
provided a good initial estimate is supplied. As we deal with relatively high SNR, this is easily realized. In 


situations where deep fades occur on the Q-band signal, a modified form of the Quinn-Fernandes 
technique is employed in which the initial frequency estimate is supplied by the K-band reference. To 
determine the Q-band frequency, a band-limited Quinn-Fernandes algorithm accurately identifies the 
beacon frequency over a reduced spectral range. In Figure 5, it can be observed that, during a rain 
event, the frequency estimator is able to track the signal as it approaches the noise floor, and recover it 
immediately once it returns. 




Figure 5. A deep rain fade event where the Q-band receiver entered frequency track mode utilizing the K-band beacon reference. 
Note that lock was maintained down to approximately -40 dB of attenuation. 
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Figure 6. Comparison of OEM tracking file-derived Doppler shift and measured frequency shift by the receivers. 


In addition, the frequency estimation also provides a means to track the relative position of the satellite 
and differentiate between signal drift and station keeping maneuvers. Plotted in Figure 6 is the deviation 
from the beacon center frequency induced by Doppler derived from the weekly provided OEM tracking 
files compared with the frequency estimator. A comparison of these two plots for 1 week show excellent 
agreement and indicate when station keeping maneuvers are taking place. Also evident throughout the 
week is a long-term mean frequency drift of approximately 15 Hz/day (K-band) and 30 Hz/day (Q-band) 
which could be a combination of system drift and/or beacon frequency drift. 


IV. Receiver Passband Calibration 

During the system checkout phase, it was observed that a diurnal drift of the signal power was evident in 
the measured data. Further investigation of this phenomenon resulted in the discovery that the filter 
passband of the IF cards exhibited some significant variation across the band. In previous campaigns, 
this effect was not noticeable, but does require calibration for the Alphasat experiment due to the non- 
trivial frequency variation induced by the Doppler shift of Alphasat’s inclined orbit. To calibrate the 
passband variation of the receivers, the antennas are pointed to zenith on a clear evening and a 
decimated version of the passband spectrum is recorded to file every second for two hours. The 
spectrum is subsequently averaged and an 8 th order polynomial is fit to the curvature of the passband, as 
shown in Figure 7. 
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Figure 7. Measurement and fitting of an 8 lh order polynomial of the passband spectrum indicating a non-trivial gain variation for (a) the 

K-DACS and (b) the Q-DACS. 


At the completion of each day, the polynomial coefficients are applied to a correction routine to generate 
calibrated data files. A comparison of the time series before and after a calibration is performed is shown 
in Figure 8, where it is observed that the calibration routine does sufficiently correct for the gain variation 
across the band. Characterization of the passband is performed nominally once a month or when major 
system changes have occurred (i.e., operational temperature changes, power outages) and the 
coefficients for calibration are updated accordingly. 
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Figure 8. Example calibration of the gain variation of the passband spectrum for (top) the K-band receiver and (bottom) the Q-band 

receiver. 


V. Data Analysis 

Preliminary analysis of the data has demonstrated excellent performance of the receiver hardware. 
Analysis of the power spectral density of a representative rain event and scintillation event are shown in 
Figures 9 and 10, respectively, and agree very well with theory. From the power spectral density plots, 
we observe that the spectrum of rain and scintillations follow extremely well with the theoretical power law 
slope of -20 dB/dec and -80/3 dB/dec, respectively [5]. In Figure 10, we observe for a clear sky day the 
effects of scintillation are evident at around 0.2 Hz and saturate at the receiver noise floor of 
approximately 1 Hz for the K-band and 2 Hz for the Q-band. A sample time series for a strong rain event 
and associated disdrometer-recorded rain rate is shown in Figure 11. 
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Figure 9. Rain event time series for 6/14/2014 (left) and the resulting power spectral density plot (right). 
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Figure 10. Scintillation event time series for 6/1 1/2014 (left) and the resulting power spectral density plot (right). 
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Figure 1 1 . Time series of a rain event for (top) the K-band receiver, (middle) the Q-band receiver, and (bottom) the disdrometer- 

recorded rain rate. 



VI. 


Conclusions 


Herein we present the design and performance of the NASA GRC-developed Alphasat beacon receiver 
terminals, which has been collecting data at the POLIMI campus since May 2014 and initiated formal 
operations on June 1, 2014. Barring system issues, the performance of the receiver has exceeded 
expectations and has been recording reliable data, which validates theoretical models of, for instance, 
rain and scintillation power spectral density slopes. We employ a modified Quinn-Fernandes frequency 
estimation routine to derive the frequency and power measurement of the beacon signals, resulting in a 
demonstrated dynamic range of approximately 38 dB at the 8-Hz sampling rate, and 45 dB at the 1-Hz 
sampling rate. The system calibration routine for passband gain variation has been described which 
adequately removes this effect. Currently, the zero-reference level calibration (not described in this 
paper) is derived from radiosonde measurements from Milano Linate Airport (approx. 7 km away). Next 
steps for this campaign involve the installation of a water vapor radiometer to more directly zero-reference 
the power measurements. Further, once the ASI site in Spino d’Adda is operational, the opportunity to 
conduct site diversity measurements at Q-band will be available, as the POLIMI site and the Spino 
d’Adda site lie approximately 20 km distance apart. 
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